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1 Introduction 



The system of partial differential equations describing a steady, irrotational and com- 
pressible fluid flow in a plane [1] is given by 

Uy - V x = 0, 

(1.1) 

(pu)x + My = 0, 

where (u, v) are the Cartesian components of the fluid velocity, and the density p is defined 
as a function of u and v. The fact that the fluid is irrotational (represented by the first 
of the two equations (11.11) ) allows us to express the velocity field in terms of a potential 
function f such that u = (p x , v = tp y . Various density functions p are admitted in the 
literature on the subject (see e.g. [21 USE]), and the most physically interesting instances 
involve bump type functions. One of the most well-investigated examples of the latter is 

p=(l + u 2 + v 2 )- 1/2 - (1-2) 

When this density function p is expressed in terms of the velocity potential <p and sub- 
stitued into the equations (jl.ip . we obtain the minimal surfaces equation in (2 + 1)- 
dimensional Minkowski space [5] 

\ I ( ^ ^ =0 (13) 

.(1 + (^) 2 + WY /2 ) X V(l + (^) 2 + WY /2 )y ' ^ ' } 

In addition, by using the Wick rotation y = it, one may transform equation (II .3p to the 
scalar Born-Infeld equation 

(l + ((fx) 2 ) <ftt - 2ip x ip t (p xt - (l - (ft) 2 ) fxx = 0, (1.4) 

which represents the simplest example of a nonlinear modification of Maxwell's electro- 
dynamics in (1 + 1) dimensions [6]. 

It is interesting to note that the Born-Infeld equation (jl.4p is compatible with the 
following hydrodynamical-type system expressed in terms of the Riemann invariants [1] 

Rt - RTRt = 0, Rt - R + R x = 0. (1.5) 
This can be shown by considering the transformation 

R ~ ± TTJ^y + TT(^- (L6) 
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Equation (11.41) is also linked with the hyperbolic Monge-Ampere equation [7] 

UxxUtt - {Uxtf + 1 = 0, (1.7) 

via the Bianchi transformations 

_ (^) 2 ~ 1 _ <fxj£t _ (^) 2 + 1 



(l-(^) 2 +(^) 2 ) 1/2 ' (l-(^) 2 +(^) 2 ) 1/2 ' (l-fo) 2 +fo.) a ) 1/2 ' 

(1.8) 

It should be noted that the system (II. 5p can be derived from the Monge-Ampere equation 
(11. 7p if we substitute 

i2 ± = — K±l). (1.9) 

We may also proceed in the opposite direction by expressing u xx and u x t in terms of R + 
and R~, and using the relations 

2R + R- R + + R~ 2 

Utt ~ R+-R- ' ^" (L10) 

It should be noted that the Monge-Ampere equation 11 1.7ft can be linearized through a 
half-Legendre transformation [8] 

u(z,y) = u(s,y) - su a (s,y), (1.11) 

where we define 

z = u s , (1.12) 
provided that u ss ^ 0. This transformation leads us directly to the linear wave equation 

Uyy-U zz = 0. (1-13) 

Finally, the system for the Chaplygin gas 

Ut = \ {U 2 - V~*) x , V t = (UV),, (1-14) 
expressed as a conservation law, can be linked to the system (jl.5p through the relations 

K k = U±± (1.15) 
and to the Monge-Ampere equation (I1.7P through the compatibility 

17=—, K = Mxx , (1.16) 
"a 

along with the relation 

w tt = c/ 2 y-y- 1 . (i.i7) 



In fact, both systems ( II. Ah and ( 11.141) can be derived, through different parametrizations, 
from the Nambu-Goto action for a string evolving in a (2 + l)-dimensional target space- 
time [I]. The Chaplygin gas (11.141) in (1 + 1) dimensions is derived from the light cone 
gauge, while the Born-Infeld equation (11.41) is obtained from a Cartesian parametriza- 
tion. The symmetry properties, subalgebra classifications and invariant solutions of these 
models have been studied extensively [U El Ell- 
in the past few years, there has also been a considerable amount of interest in the 
study of theories involving odd (fermionic) Grassmann variables and super algebras. Such 
systems are interesting because ordinary matter generally consists of fermions wheras 
bosons are only concerned with the interactions. One relatively recent line of inquiry in 
this area of research has involved the construction of supersymmetric extensions of existing 
classical and quantum systems. This approach, first used in the context of particle physics, 
was successfully employed to super symmetrize theories involving classical and quantum 
fields jH[Lll[T2]. Application of the techniques to fluid dynamics began with the study of 
simple Euler-type systems [T3| [T4"l [T5] , followed by extensive work on a supersymmetric 
version of the Korteweg-de Vries equation [TH [TTJ, [18] . More recently, the Chaplygin gas 
was supersymmetrized in both (1 + 1) and (2 + 1) dimensions by R. Jackiw, Y. Bergner 
and A. P. Polychronakos [H [T9], [20] through the addition of fermionic-valued fields to 
the classical theory in the bosonic field 9. In those models, the velocity was no longer 
irrotational but was expressed in terms of both the bosonic and fermionic fields 

v = W-|^VV>. (1.18) 

Their approach involving a Lagrangian formulation of the fluid dynamics equations has 
been applied, among other areas, to an N = 1 supersymmetric extension of polytropic 
gas dynamics [21], a covariant and supersymmetric theory of relativistic hydrodynamics 
in four-dimensional Minkowski space [221 123], and a Kaluza-Klein model of a relativistic 
fluid [21]. 

It must be said that the physical interpretation of the systems resulting from the 
supersymmetric extensions remains an intriguing and still open question. However, a 
number of remarkable, physically meaningful solutions obtained in both particle physics 
[TT1 [25] and fluid dynamics [26J attest to the potential of this methodology and motivate 
further attempts at its application. In this paper, we apply it to a Gaussian fluid flow 
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which, to our knowledge, has not been supersymmetrized before. 

The objective of this paper is to investigate the equations of a steady, irrotational and 
compressible fluid flow (11.11) involving the Gaussian density function 

p = e- u2 - v \ (1.19) 

This type of Gaussian fluid plays an essential role in many areas of physics, including fluid 
dynamics [3l [27], plasma physics [28] and astrophysics [3]. When the density function 
(I1.19P is substitued into the equations (11.11) and the velocity components expressed in 
terms of the potential <p, we obtain the equation 

(l - (f x ) 2 ) (p xx - 2(p x cp y (f xy + (l - {Vy) 2 ) Vyy = 0- (1-20) 

Through the use of a method similar to that used for the Born-Infeld equation (11.41) . we 
will formulate a supersymmetric generalization of equation (ll.20p . Due to its interesting 
symmetry properties and solutions, we will also consider a slightly modified version of 
equation (11.201) found in statics of liquids [29J . 

(l + [ip x ) 2 ) if xx - 2(p x ip y (f xy + (l + (y?y) 2 ) <Pyy = 0. (1.21) 

For the purpose of considering both equations (ll.20p and (ll.2ip . let us introduce the 
parameter e = ±1. The values e = 1 and e = —1 will correspond to equations (ll.2Up and 
(ll.2ip respectively. Thus, the two equations may be represented in the form 

(l - e{ip x ) 2 ) ip xx - 2ip x ipy(p xy + (l - e{ip y ) 2 ) ifyy = 0. (1.22) 

The goal is to construct supersymmetric extensions of both the Gaussian fluid flow equa- 
tion (ll.20p and its modified version (11.211) . to study their symmetry properties and to 
obtain new classes of invariant solutions of the extended models. The approach used 
in this paper is different from that used by Jackiw, Bergner and Polychronakos for the 
Chaplygin gas, and is based on the symmetry reduction method adapted to Grassmann- 
valued partial differential equations. It has been used in the past for the supersymmetric 
Korteweg-de Vries equation [T5] and, more recently, by one of the authors for the super- 
symmetric Born-Infeld scalar model [50] . 

Our paper is organized as follows. In section 2, we identify the symmetries, subalgebra 
classifications and invariant solutions of the classical Gaussian fluid flow equation (11.221) . 
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In section 3, we describe the most general form of the supersymmetric extension in terms of 
a superspace formalism. In section 4, we examine the symmetry properties of the extended 
system and compare them to those of the supersymmetric scalar Born-Infeld equation 
constructed previously and describe the subalgebra classification structure. In section 5, 
we use the symmetry reduction method to obtain invariant solutions, including elementary 
solutions (algebraic with one and two simple poles, trigonometric and hyperbolic) and 
doubly periodic solutions which can be expressed in terms of elliptic functions. Finally, 
section 6 contains a summary of our results. 



2 Symmetry properties of the classical Gaussian fluid flow 



Before we construct our supersymmetric extension, we first examine the Lie point sym- 
metries and invariant solutions of the classical irrotational Gaussian fluid flow (jl.22p . 



2.1 The case where e — 1 

We consider first of all the case where e = 1. The Lie symmetry algebra of equation 
(jl.20p is spanned by the following five generators 

S(i) = xd x + yd y + y?<9 v , M = -yd x + xd y , 

(2.1) 

T\ = d x , T 2 = d y , T3 = d v , 

which represent a dilation, a rotation and three translations in each of the three vari- 
ables. The commutation relations of the algebra spanned by the generators (12.11) are 
summarized in table 1. 

Table 1: Commutation table for the Lie symmetry algebra Qn\ generated by the vector fields 





S(i) 


M 


Ti 


T 2 


T 3 


S(i) 








-Ti 


-T 2 


-T 3 


M 








-T 2 


Tt 





Ti 


T X 


T 2 











T 2 


T 2 


-Ti 











T 3 


T-i 
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The Lie algebra Q(\) can be decomposed as the semi-direct sum 

= {S(i),M} 35 {Ti,T2,T 3 }. (2.2) 

We now classify the one-dimensional subalgebras of the Lie algebra Qm using the clas- 
sification methods as described in [31 J . That is, we construct a list of representatives of 
conjugacy classes of subalgebras in such a way that each one-dimensional subalgebra is 
conjugate to one and only one element of the list. We focus on one- dimensional subalge- 
bras with orbits of codimension 1 because we are only interested in invariant (as opposed to 
partially invariant, conditionally invariant or generic) solutions, so that one-dimensional 
reductions will lead to ordinary differential equations. For each of the elements of the 
following list, we first identify the representative subalgebra of the conjugacy class. Next, 
we describe the invariants of the subalgebra and (except in the case of Li) the change 
of variable. Finally, using the symmetry reduction method as described in [32l 133] . we 
substitute the change of variable into the original equation f ll.20p in order to obtain a 
reduced ordinary differential equation in the symmetry variable £. 

L\ = {T 3 }, Invariants: x, y; 

L 2 = {Ti}, Invariants: £ = y, (p, Change of variable: (p = <p(y), 

(1 - (ip y ) 2 ) ip yy = 0; (2.3) 

L 3 = {T 3 + mTi} m ^ , Invariants: f = y, <p - —x, 

m 

Change of variable: 02 = F(y) H x, 

m 

(1 - (F y ) 2 ) F yy = 0; (2.4) 

x 

Li = {S(i)\, Invariants: £ = — , — , Change of variable: ip = yF(£), 

y y 

[(1 + e) - (1 + a 2 (^) 2 + 2£(i + e)FF^ - ep 2 ] % = o ; (2.5) 

L 5 = {M}, Invariants: £ = x 2 + y 2 , cp, Change of variable: (p — <p(£), 

<Pt + &k - 2^ n ) 3 - 4^ 2 (^)V« = 0; (2.6) 



L 6 = {Sqj + aMj^o, Invariants: £ = arctan ^— J — a In {^J x 2 + y 2 ^j , 



y/x 2 + y 2 



Change of variable: (p = \/ x 2 + y 2 F(£), 

-2a(l + a 2 )FF^ + (1 + a 2 ) 2 (F^ 2 F^ + a 2 F 2 F^ - (1 + a 2 )F^ (2.7) 
-a(l + a 2 )(F^ 3 + (1 + 2a 2 )F(F^) 2 - aF 2 F^ + 2aF^ - F = 0; 



/ y 

Lj = {M + [iTsj^i, Invariants: £ = x 2 + y 2 , (p — /zarcsin 



\/ x 2 + y 2 



Change of variable: ip = F(£) + /iarcsin 



y 



V x2 + y 2 



F c + 2^ - 4^(F c ) d + 2^F^ - 8r(F € ) 2 % = 0. (2.8) 

In the case of the subalgebra L%, we obtain no reduced equation and therefore no invariant 
solution. This corresponds to the fact that there is no function (p(x,y) which is invariant 
under the transformation ip — > ip> + K, where K is a constant. Solutions invariant under 
the subalgebras L 2 and L 3 are linear polynomials in x and y. 

For subalgebra L 4 , we obtain the following three solutions for equation (12.51) 

F(0 = C£ + C 2 , (2.9) 

which is a simple linear function, and the two expressions 

= v 7 ! + (± arctan f + Ci) , (2.10) 

which correspond to a kink-type fluid density (a condensation wave). For the subalgebra 
L e , the solution of equation (12.61) is 

1/2 



0(0 = ±- J I y 7 I + C 2 , (2.11) 

where the Lambert function, y = C(x), is the solution to the equation 

ye y = x. (2.12) 

Although, for each x, there are infinitely many values which satisfy (I2.12p . it should be 
noted that exactly one of the branches is analytic at zero. The multi-valuedness of this 
solution exhibits ergodic behavior [34J. 



Finally, for subalgebras L 6 and L 7 , equations ( 12 .TP and (12.8!) do not possess the Painleve 
property. 
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2.2 The case where e = 



-1 



Considering now the case where e = — 1, we see that the Lie symmetry algebra G(-i) of 
equation (11.211) is spanned by the following four generators 

= xd x + yd y + ipdp, t x = d x , * 2 = d y , * 3 = d v , (2.13) 

which differs from the case e — 1 in the sense that there is no rotation present. It should be 
noted that this algebra is solvable. A classification of the one-dimensional subalgebras of 
G(-i) leads us to conclude that each subalgebra of the form {0*1 + 6*2 + 0*3} (where a, 6 and 
c are real constants) is in its own separate conjugacy class, whereas every other subalgebra 
is conjugate to {S(_i)}. Using the classical symmetry reduction method, one can obtain 
the following invariant solutions of ( 11.21 j> . A solution invariant under the action of the 
dilation SV_i) must be of the form (p = yF(£), where £ = x/y is the symmetry variable. 
In this case, equation (jl.2ip reduces to the ordinary differential equation 

((1 + i 2 ) + (1 - £ 2 ) 2 (^) 2 + 2f (1 - £ 2 )FF ? + ?F 2 ) F u = 0. (2.14) 

Due to the factorization of (|2.14|) . one obtains directly a solution which is linear in £ 
and corresponds to the (linear) solution u(x, y) = C\X + C 2 y, where C\ and C 2 are real 
constants. In addition, one obtains the following two solutions expressed in terms of 
elliptic functions 



F(0 =C lV ^Ty^+l 



(2.15) 



where 



™ = f ji-M-ec?*' (2 ' 16) 







is the incomplete elliptic integral of the first kind (the inverse of the Jacobi function sn) 
and 



a 



£{z,k)= —, da (2.17) 







vT — ( 1 



is the incomplete elliptic integral of the second kind [35]. Solutions invariant under a 
subalgebra of the form {a*i + 6*2 + 0*3} all consist of linear polynomials in x and y. 



9 



3 Supersymmetric extension 



In analogy with the case of the supersymmetric Born-Infeld scalar equation, we construct 
a supersymmetric extension of equations ( 11. 22ft through a superspace Grassmannian for- 
malism. More specifically, we supplement the set of bosonic (even) independent variables 
{(x,y)} with a fermionic (odd) Grassmann parameter 9 and replace the bosonic function 
<p(x, y) with the superfield <&(x, y, 6) defined as 

$(x,y,e) = il;(x,y) + 6(p(x,y), (3.1) 

where ip(x,y) is a fermionic field. We construct our extension in such a way that it is 
invariant under the supersymmetry transformation 

x^x-yd, e^e + v, (3-2) 

where rj is a constant fermionic parameter. This transformation is generated by the 
infinitesimal supersymmetry operator 

H = d e -6d x . (3.3) 

In order to make our superfield theory manifestly invariant under the action of H, we 
write the equation in terms of the covariant derivative operator 

D = d e + 9d x , (3.4) 

which possesses the property that it anticommutes with the operator H. That is, 

{H, D} = ED + DH = 0. (3.5) 

The most general form of a supersymmetric extension of (ll.22p is given by the fol- 
lowing equation in terms of the superfield $ and various orders and combinations of the 
derivatives D and d y 

D% - Ea(D 2 ^)(D :i ^)(D 5 ^) - e(1 - a)(D 3 $) 2 (D 4 <I)) - 2b{D 2 <$>)(D§) y (D*<$>) y 

- 2c(D^)%(D^) y - 2(1 - b - c)(D 3 $)(D$) y (D 2 $) y + & yy (3.6) 

- ed<fry(D&)y(D&)yy - 6 (l - d) {{D*)yf ^ = 0, 

where a, b, c and d are arbitrary real parameters. Here, the subscript y indicates differen- 
tiation with respect to y. This equation can be decomposed into the following two partial 

10 



differential equations for the fields (p(x,y) and i[)(x,y) corresponding respectively to the 
coefficient of the fermionic variable 9 and that of the remaining terms 

Vxx - £(¥x) 2 ¥xx + £aip xx ip x ip xx + eaip x ip x ^ xxx - 2(p x ip y ip xy + 2b(p xy tjj x il) X y + 2b(p y ip x ijj XX y 

+ 20p x y1p y 1p XX + 20p x ll)yll) XX y ~ 2 (fy ljj xx lj) x y + 2l)(Pylj) xx lj) X y + 2(Xf >y1p ' XX 1p ' X y + ^ yy 
- e(iPy) 2 (fyy + B 'tfyj <yl/) > ylf) ' X yy + Sdif yyl^ yljj X y ~ 26 if y^ X y^ yy + 26 (If y^ X y^ yy = 0, 

(3.7) 

and 

%j) xx - eaif x if xx ip x - e(ip x ) 2 i! xx + ea(tp x ) 2 ip xx - 2bip y if xy i/) x - 2ctp x tp xy i/;y - 2ip x ifyif) X y 

+ 2bip x iPy4> X y + 2Cif X iPyi) X y + lf)yy ~ Sdlfylfyy^y ~ S(<f < y)' ''if, ) yy + S d(if ' y)' ' l[) ' yy = 0. 

(3.8) 

In the limit where ip — > 0, equation (13.71) reduces to the standard bosonic equation (11.221) . 
Equation (13.81) is a new equation involving the fermionic field if). 



4 Lie symmetry superalgebra 

For the system of partial differential equations defined by the system (13. 7p and (13.81) . the 
following symmetries are present for all values of the parameters a, b, c and d. We have 
the following dilation in independent and dependent variables 

S = xd x +yd y + <pd v + lil>d 1 i„ (4.1) 

along with the following four translations 

Pi = d x , P 2 = d y , Z = d v , F = <V (4.2) 

In addition, in the case where a = and b = 0, we have the following additional fermionic 
symmetry 

Qi = a% (4.3) 
Similarily, in the case where c = and d = 0, we have the additional generator 

Q 2 = yd*. (4.4) 

Finally, it should be noted that the generators H and D given in equations (13.31) and 
(13 .4p . which involve the independent variables x and 9 in the superspace, correspond to 
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the following generalized symmetries of the coordinate space X X U = {(x, y, <p, ip)} 

H = -ipxdp + ipd$, (4.5) 

and 

D = i) x d v + (pd,p. (4.6) 

It is interesting to note that the symmetry generators found for our supersymmetric 
Gaussian fluid flow equations (13.71) and (13.81) are analogous to those found for the su- 
persymmetric scalar Born-Infeld extension [3D]. The only instance for which all seven 
symmetries (14. ip . (14.21) . (14.31) and (14.41) are present is the case where a = 0, 6 = 0, c = 
and d = 0. From here onwards, we shall focus our attentions on this specific case. The 
equations (13. 6p . (13.71) and (13.81) become 

D% - e(D^) 2 (D%) - 2(D^)(D<b) y (D 2 <b) y + % y - e {{D<£) y f % y = 0, (4.7) 

VXX - £(<Px) 2 <PxX - 2ip x (p y ip xy - 2(p y ljj xx ljj xy + Ifiyy ~ 6 (if ' y)' 'if ' yy ~ 2Slf ylj) x ylj) yy = 0, (4.8) 

and 

tP XX - e{lf X ) 2 l\) XX - 2lf X lfytjj X y + tjjyy ~ E (if > y)^ ^ 'yy = 0, (4.9) 

respectively. The commutation and anticommutation relations of the superalgebra L 
formed by the generators (14.11) . (14.21) . (14.31) and (14. 4p are summarized in table 2. 

Table 2: Commutation table for the Lie symmetry superalgebra in the case where a, b, c, d = 





s 


Pi 


p 2 


z 


Y 


Qi 


Q 2 


s 





-Pi 


-P2 


-z 


2 1 


-5Q1 


-5Q2 


Pi 


Pi 














Y 





p 2 


P2 

















Y 


z 


z 




















Y 


2 1 




















Qi 




—Y 

















Q 2 


5Q2 





-Y 















4.1 One-dimensional subalgebras 

In this section, we describe the classification of the one- dimensional subalgebras of the 
Lie superalgebra L obtained in section 3. In other words, we proceed to construct a list of 
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representatives of conjugacy classes of subalgebras in such a way that each one-dimensional 
subalgebra is conjugate to one and only one element of the list. This classification has 
already been performed for the Born-Infeld superalgebra [30], and so we describe the 
results. We focus on one-dimensional subalgebras with orbits of codimension 1 because we 
are only interested in invariant (as opposed to partially invariant, conditionally invariant 
or generic) solutions, so that one-dimensional reductions will lead to ordinary differential 
equations. The analysis led to the following results. 

The one-dimensional splitting subalgebras are 

L x = {S}, L 2 = {P 1 }, L 3 = {P 2 }, U tm = {Pi + m? 2 }^ , 

L 5 = {Z}, L 6jTn = {Z + mPi} m ^ , L 7jm = {Z + mP 2 } m ^o, 

L 8 ,rn,n = {Z + mP x + nP 2 } m>r ^ , L 9 = {Y}, Lio = {Qi}, L n = {Q 2 }, 

L\2,k = \Qi + kQ 2 }k^o, 

(4.10) 

and the non-splitting subalgebras are 

£2,rn,w = { p i + ViQi + V2Q2}, ^3,rn,m = { p 2 + V1Q1 + V2Q2}, 

£4,m,rn,V2 = {Pi + mP 2 + VlQl + ??2<22}m^0, £>5,m,V2 = { Z + VlQl + V2Q2} , 

£>Q,m,7n,m = {Z + mP\ + V1Q1 + Wo, I 4 - 11 ) 

,771,771,172 

{Z + mP 2 + rjiQi + Wo, 

£8,771,71,771,772 = { Z + m -P + nP 2 + VlQl + ^^jm^O, 

where m, n and k are constant bosonic parameters, and r]i and r] 2 are fermionic constants. 
Throughout this paper, we use the convention that underlined constants are fermionic. It 
is important to note that each non-splitting subalgebra Ci jVliV2 is an element of the same 
conjugacy class as the subalgebra C^k^Kt^, where K is a positive constant. 



5 Symmetry group reductions and solutions 

In this section, we make use of the classical symmetry reduction method to determine 
the invariants and reduced equations corresponding to each of the one-dimensional sub- 
algebras in equations (14.101) and (14.111) . In each case, we obtain at least one invariant 
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expressed solely in terms of the independent variables x and y, which we call the symme- 
try variable £. In most other cases, the remaining two invariants involve the dependent 
functions (p and ip, and they are identified as the functions F(£) and A(£) respectively. 
Expressing the fields tp and ip in terms of the invariants, we substitute the derivatives 
into the supersymmetric system (14 .8p and (14. 9 p in order to obtain the reduced a system 
of ordinary differential equations. The results are presented in tables 3 and 4. Where the 
Painleve property is satisfied for these ODEs, we give an exact, analytic solution of the 
system of reduced equations, along with the corresponding solution of the supersymmetric 
system ( 14. 10D and ( 14. lip . It should be noted that this reduction procedure cannot be em- 
ployed for subalgebras whose symmetry generators do not involve derivatives with respect 
to independent variables, nor for the generalized symmetries (14 .5p and (14. 6p . Indeed, such 
transformations involve shifts only in the dependent variables, which cannot under any 
circumstances leave a solution invariant. 

5.1 Splitting subalgebras and their reductions 

We begin our analysis by considering the splitting subalgebras. In the case of subalgebra 
Li, we obtain two equations involving the functions F and A, coupled in an involved way. 
In the specific case where the condition 

(f eAA 5 - | eZAAx + (1 + ^ 2 )A C A 5C ) (F - £F ? ) = (5.1) 

holds, we can decouple the first equation which then becomes an ordinary differential 
equation in F 

((1 + i 2 ) - e{\ + eefW 2 + 2f (1 + ee)FFt - ^ 2 F 2 ) F^ = 0. (5.2) 

In the case where e — 1, there are three functionally independent solutions of equation 
(15. 2p . The first two are expressed in terms of radicals and inverse trigonometric functions 

F{0 = Vl + f 2 (±arctanf + Ci) (5.3) 

Here, the function A can be determined from the second reduced equation in combination 
with the condition (15. ip . and is given by the expression 

A(f) =Ei (1 + £ 2 ) 3/4 (±arctan£ + Ci) 3/4 exp ^ - (±|Ci + § arctan£) arctan^j . (5.4) 
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Table 3: Invariants and change for variable for each splitting and non-splitting subalgebra 



Subalgebra 


Invariants 


relations and change of variable 


Li — W\ 




10 — vF(f) ib — ij 3 / 2 A(£) 


L 2 = {Pi} 




<p = v(y),ib = ^(y) 


^3 — X* 2/ 


op I r\ 1/1 


ip — <p{X),ip — tp{X) 


£4,771 = {Pi + mP 2 } m ^0 


i = V- mx,ip,i/j 


f = v(0.^ = ^(0 


r ; , — fz\ 


T* 7/ 7/1 




L 6im = {Z + mPi} m ^ 




^=^(^(^+^^ = ^(2/) 


^ ( ,771 1 Z J TTl^U 


X. 7/ — 77109. 1/) 


03 = — (Fix) + v) , ib = ib(x) 


-^8, 771,71 — 

{Z + mPi + n-P^i™ 


t = x- ?%y,<p- %,ib 


<P = F(0 + ±y,*b = ib(0 


L 9 = {Y} 


x,y,ip 


N/A 


Lio = {Qi} 


x, y, if 


N/A 


Ln = {Q 2 } 


x,y,ip 


N/A 


Li2,k = {Qi + kQ 2 }k^o 


x,y,ip 


N/A 


^2,771,7)2 — {-^1 
+V1Q1 + V2Q2] 


y,ip,ip - \mx 2 - mxy 


cp = <p(y),ib = A(y) + ^rji x ' 2 + "mxy 


£3,771,772 = {-^2 

+ %(32} 


x,ip,ip- rjixy - fyvy 2 


tp = <p(x),ib = A(x) + rjixy + ^my 2 


£4,771,771,772 = {A + "1-P2 
+77lQl + ?72<22 Wo 


i = y- mx,ip, 
ib — \r\ix 2 + \r} 2 mx 2 — i] 2 xy 


f = ¥>(£)> 
V' = A(^) + ^r/ix 2 — \r} 2 mx 2 + r72xy 


£5,7)1,7)2 = {-^ 
+?7lQl + ??2<22} 


x,y,tp- {mx + my) v 


AT/A 


£6,711,771,7)2 = + mPi 
+ ?72<22 Wo 


y,<p-±x,il>- ^mx 2 - ^mxy 


v = ^(y) + fb> 


£7,771,771,7)2 = + "2-^2 
+*7lQl + <72<22 Wo 


XiV-^y^- imxy - ^ W 


V' = A(x) + ±vixy + 2^%y 2 


£8, 771, 77, 7)1,7)2 

{Z-r-mP7TnP 2 

+77lQl + r?2Q2}m,n^0 


Z = x- ^y,<p-±y, 
ip - i\m_xy + ^m 2 - ^my 2 
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Table 4: Reduced Equations corresponding to each splitting and non-splitting subalgebra 



Subalgebra 


Reduced equations 


Li = {D} 


(1 + a% - e(i + eey{F^F^ + 2^1 + ee^F^ 

-e(?F 2 Fc i + |e(F - £F,)AA e - §ef (F - £F )AA ee 
+ (1+^ 2 )(F-^)A ? A K = 0, 

(1 + £ 2 )A ?? - £ (1 + e£ 2 ) 2 (P?) 2 A ff + 2C(1 + ee)FF^ 
-e£ 2 P 2 A K + £(1 + ^ 2 )(F f ) 2 A c - £(2£ 2 + e)FF^ 
+ e £F 2 A f - £A 5 - f ee 2 (P?) 2 A + §e£FF c A - | eF 2 A + § A = 


L 2 = {Pi} 


- £(fy) 2 fyy = 0, 

Ipyy ~ e(ip y ) 2 1p yy = 




- s((p x ) 2 ipxx = 


La m — {-Pi + ^P^Im^O 


(m 2 + 1W - £(m 2 + £) 2 (w f ) 2 o9 ff = 0, 
(m 2 + - e(m 2 + e) 2 (^) 2 % = 


L 6 ,m = {Z + mPi} m ^ 


^ Pyy ^(Py) Pyy — 
™ 2 ^ra - ^(Pj/) 2 ^ = 


L>7,m = {Z + 771,^2 }m^0 


e{Fx) F xx = 0, 
m 2 ^ - e(F x ) 2 ip xx = 


L,o = \Z + TTlPl + ni^oTm n^n 

0,7/1,71 L 1 i 1 Z J Hijfl^U 


f_ e (l + e ™i) 2 (F0 2 + (2^ + 2e^t)F f + (1 + - eK)) F ee = 0, 
(-e(l + e£) W + (2$ + 2 £ ^)F ? + (1 + £ - e £)) = 


•£2,771,772 = {-Pi 


-2eipyr]2Ayy - 2^r?ir? 2 + - e((p y y 2 <p yy = 0, 
A ra - e(^) 2 A ra + m = 


£3,771,772 = { p z 
+mQi + 1I2Q2} 


(l - e^) 2 ) Vxx = 0, 
A ra - e{ip x ) 2 k xx + 772 = 


£4,771,771,772 = {-Pi + 


(1 + m 2 )ip^ - e(l + em 2 ) 2 ^) 2 ^ 
+2<^ ((-£772 + mr/i)A §5 + 772771) = 0, 
(1 + m 2 )A K - £(1 + £m 2 ) 2 (<^) 2 A^ 
+m(2r]2 - £m7/i + em 2 r]2){(p^) 2 + (771 - 771772) = 


£6,771,771 ,772 = {-^ + 777,Pl 


Fyy - e(Fy) 2 Fyy ~ ^Ti^F,, - E-mFykyy = 0, 

A ra - e(Pj / ) 2 A ra - ^2-772^2/ - £^(1 - £m 2 )7/i = 


£7,771,771,772 = {Z + mP2 
+V1Q1 +mQ2}mjtO 


Prx - e(F x ) 2 F xx + ^2-771 A ra - e^mm = 0, 
A^ - e(F x ) 2 A xx - ^mFx - £^r(l - etti 2 )^ = 


£8,771,71,771,772 ={Z + mPi + nP 2 
+ %Q2}m,n^0 


(-£(1 + £^) 2 (F f ) 2 + (2$ + 2£^)F $ + (1 + ^ - £^)) % 

+^ (m " (1 " mF^A^ -e£(l- mF^ mm = 0, 
(-£(1 + £^) 2 (F C ) 2 + (2^ + 2£^)F 5 + (1 + £ - % 
+ ("^^ " + (F,) 2 + £^(F ? ) 2 - ^ + e%) 771 
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Here, condition ( 15. ip is identically satisfied due to the fact that the constant E\ is 
fermionic. The solution (15.31) corresponds to a kink-type fluid density 

p = e -(l+(-ctan f+ C 1 ) 2 ) ; (55) 

where the velocity components are given by 



farctanf +d) + ^= (5.6) 



\Jx 2 +y 2 V v J ^x 2 +y 2 

v = . I ( arctan 5 + - , x . (5.7) 

-fi+f— — +c ) 2N ) 

This is a kink-type density solution, which varies from the value e \ ^ 2 ^ / to the 

_(l + (n +c \ 2 \ 

value e V ^ 2 1 > and depends only on the polar angle of the position on the plane. 
This solution is real, asymptotic and discontinuous on the x-axis, and represents a con- 
densation wave. 

The third solution of equation (15. 2p is the linear function 

F(0 = C^ + C 2 , (5.8) 
which gives us the following solution for A 

A (0 == 3/2 - (5.9) 

((1 - CU + C X C 2 + y/Cl + CI - l) 

In the case where e — — 1, there are again three independent solutions of equation (15.21) . 
The first two are given in terms of elliptic functions 

F{0 =c iy /tr^iy/i+i± 

where C\ G R, and fc) and fc) are the incomplete elliptic integrals of the first and 
second kinds as defined in equations (12. 161) and (12. 17ft respectively. Solution 115. 10p is a 
complex- valued function, but for large positive values of C\, the imaginary part provides 
a good approximation for a localized static bump solution. The function A is given by 
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(5.11) 



the expression 

A(0 exp f ? J (c, i) - % 0) (2^° + 6£ 8 + 4£ 6 - 4£ 4 - 6£ 2 - 2) 

+ V z ^ r e^ + a 7/2 v^a/^+Vk - m - 1) ■ 
(l+c^+^+^e.o-^e.o) 3 )) / 
(te - i)(e + 1) ( (£(£, o - o) a/ 3 ! 3 ^ vTh (i + e 2 ) 4 

+c 1 (i+a 9/ V(-e-i)(^-i) 

+ x^rq^i + e 2 ) 7 / 2 v / ^T v / e+Tv / (-^ - W - 1)))^ • 

The third solution of equation ( 15. 2ft is the linear function 

F(fl = Ctf + C 2) (5.12) 
which gives us the following equation for A 

(1 + C 2 X + 2C,C 2 i + e + Cf£ 2 )A« + {-C 2 2 i - i - dC 2 )A e + |(1 + C|)A = 0. (5.13) 
Solving equation (15.131) . we obtain the hyperbolic solution 



= c ^ sinh (/ 2^i + cj g (o * 

+ Cl vw> cosh f / s_vi±s±^^^mi±^A , 



(5.14) 



where we have defined 

<7(0 = 3C 2 C 2 2 - C 2 - C? - 1 + 6dC 2 e + QClC^ + 6C 2 £ 2 + 3C 2 4 £ 2 + 3£ 2 . (5.15) 

For subalgebra L 2 , we obtain the following two solutions. First, a trivial linear solution 
in y for both if and ^, and second a solution which is linear for ip (<p(y) = ±iy + C 2 ) but 
where ip is an arbitrary function of y. In particular, this includes the case where ip is a 
solitonic, bump, kink or multiple wave solution. From subalgebra L3, we obtain solutions 
similar to those for L 2 except that the argument y is replaced by x. 

For subalgebra L4, we should distinguish two separate cases. In the case where e — — 1 
and m = ±1, we obtain the linear travelling wave 

<p(x,y) = C 1 (y±x) + C 2} ^(x, y) = Ki(y ± x) + K2, (5.16) 

18 



where C\ and C 2 are bosonic constants and K\ and K 2 are fermionic constants. For all 
other cases of e and m, we obtain two separate families of solutions, one of which is linear 

(p(x,y) = Ci(y-mx) + C 2 , i/>(x, y) = Ki(y - mx) + K2, (5.17) 

and the other, where tp(x,y) is given by the specific expression 

/ m 2 + 1 \ 1/2 

<p(x,y) = ± \£ ^ m2 + £ y ) (y~mx) + C 2 , (5.18) 

while the fermionic field ip is an arbitrary function of the single quantity y — mx. This 
allows us to consider a wide range of interesting physical phenomena, including travelling 
waves, center waves, bumps, kinks and multiple waves. In particular, we can choose 
fermionic solitary wave solutions, which are particularly interesting since they have been 
extensively studied and could give us information about the fermionic medium under 
consideration. 

From subalgebras L$ and L7, we obtain linear solutions for (p similar to those discussed 
above for L 2 and L 3 above, respectively. We still obtain solutions for the fermionic field 
ip which are arbitrary functions of y and x respectively. For subalgebra L 8 , we get yet 
another linear solution similar to that of L 4 , but also the specific solution 

<p(x,y) = \ ; 2 J -(x--y)+C 2 , $ = $(x--y), 5.19 

m 2 + en 2 \ n / \ n / 

in terms of an arbitrary function for ■0. This represents a propagation wave in x and 
y. Because of this freedom in ip, we can consider, for example, bounded elementary, 
trigonometric, periodic and doubly periodic solutions, as well as Painleve transcendents. 
Consequently, different physical phenomena can be considered. 

5.2 Non-splitting subalgebras and their reductions 

Let us turn now to the non-splitting subalgebras. Combining the two reduced equations 
for subalgebra C 2 , we obtain the following second order ordinary differential equation for 
ip only 

(1 - e(ip y ) 2 ) 2 Vy y - ( y9 y ) 3 + 4^,1^1^ = 0, (5.20) 

where 

{1, if a = b 
(5.21) 
0, if a ^ b 

19 



is the usual Kronecker delta function. Equation (I5.20j) can be converted through the 
substitution to = ip y to a first order equation 

(f - s(uj) 2 ) 2 uj y - 2TJITJ2UJ 3 + 4^,1^771772 = 0. (5.22) 

In the case where e — — 1, equation (15.221) can be integrated to give the transcendental 
relation for u 

4u 2 \nu- Arjjjfr^y + - 1 - Adrj^u 2 = 0. (5.23) 

The subalgebra £3 leads us to the quadratic solution 

1 

2(1 -eC\ 

In addition, in the case where r] 2 = 0, we also have 



<p(x) = G x x + C 2 , i[)(x) = - — -p^i^x 2 + + K2 + rjixy + \rny 2 . (5.24) 



(f(x) = ±y/ex + C 2 , ip(x) = A(x) + rjixy, (5.25) 

where A is an arbitrary function of x. In similarity with the previous cases, we can adjust 
A to various solitonic, compact suport functions. 

The two reduced equations for subalgebra £ 4 can be combined into the following second 
order ordinary differential equation for (p 

e(l + m 2 )V« - 2(1 + m 2 )(l + em 2 ) 2 (<pe) 2 (p# + e(l + em 2 ) 4 (p e )V« 

(5.26) 

+ 2(1 - m 2 )r/i?72(v9 ? ) 3 - Mj^tp fflrfr = 0, 
which can be reduced to the first order equation for u 
e{\ + m 2 ) 2 u, - 2(1 + m 2 )(l + em 2 ) 2 uj 2 uj, + e(l + em 2 ) Vuv + 2(1 - m 2 )^^ 3 

(5.27) 

- 46 £jl urjiT]2 = 0. 

This equation does not have the Painleve property, implying that there are additional 
singularities besides fixed poles. 

For the case of subalgebra Cq, the reduced equations from table 3 combine to give the 
following equation for F 

(1 - e(Fy) 2 ) 2 F yy + e^rjMFy) 3 + -^(1 - 2m%, 1 )r m F y = 0. (5.28) 

This equation can be transformed to the first order ordinary differential equation for uo 

2 2 2 

(l-euj 2 ) uj y + ri^u 3 + — (1 - 2m 2 5 eyl )r] l 7] 2 uj = 0. (5.29) 
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For the case of subalgebra £7, combination of the reduced equations gives us the 
following equation for F 



(1 - e(F x ) 2 ) 2 F xx + ^mrn (m 2 (F x ) 2 + e - 2m 2 8 e>1 ) = 0. (5.30) 



This equation corresponds to the first order ordinary differential equation 

2 

(1 - eu 2 ) 2 u x + —771772 (mV + e - 2m 2 5 eA ) = 0. (5.31) 



For subalgebra £ 8 , the two reduced equations give us the equation 

(m 2 + en 2 ) 4 (F^) 4 - Aem(n 2 + em 2 ) 3 (^) 3 

+ (Qm 2 (m 2 + en 2 ) 2 - 2en 2 (n 4 + em 4 )(n 2 + em 2 ) - 4m 2 n 4 (n 2 + m 2 )5 £ ^ (F ? ) 2 
+ Am{n 2 {n i + em 4 ) — em 2 {n 2 + em 2 ) + 2m 2 n 4 5 £ ^F^ 

+ (m\l-en 2 )-en 4 f^ F^ (5 - 32) 

+ ^ - 2emn 3 (m 2 + en 2 )(F ( f + 2n 3 (n 2 + 3em 2 )(F ( ) 2 

+ (4mn 5 4,i - 6emn 3 )F^ + 2er? - 4n 5 5 £ ^ 771772 = 0. 
Through the substitution u = ip y , equation (15.321) can be reduced to a first order equation. 

6 Conclusions 

In this paper, we have formulated supersymmetric generalizations for the Gaussian irro- 
tational fluid flow equation (11.201) and its modified version (11.211) . It is interesting and 
significant to note that the symmetry superalgebras for our supersymmetrized Gaussian 
fluid model equation and for its modified version were isomorphic to each other, as well 
as to that of the supersymmetric Born-Infeld scalar equation (11.41) as investigated in [30] . 
This is in contrast to the classical case, where the symmetries for equation (jl.20p included 
a rotation which was not present in (ll.2ip . The three translations in x, y and present 
in the classical case are also present in the supersymmetric case, while the dilation in 
x, y and <fi is modified to a dilation which also includes the fermionic field ip. However, 
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the rotation transformation present in the classical fluid equation ( 11.201) is not preserved 
in the corresponding supersymmetric equations (13 .7p and (13.81) for e — 1. Conversely, a 
new translation in ip and two more fermionic transformations (for specific values of the 
parameters a, b, c and d) are present in the supersymmetric case. 

Through the use of the symmetry reduction method, invariant solutions of the classical 
and supersymmetric Gaussian fluid were systematically constructed. Solutions of the 
classical Gaussian flow equation fll.20j) include a kink-type fluid density (I2.10p representing 
a condensation wave and a solution (12.111) expressed in terms of the Lambert function. The 
kink-type expression (I2.10p is also present as a solution of the supersymmetric equations 
(14.81) and (14.9[) with e — 1, where it represents the bosonic part of the solution (15. 31) and 
(15 .4p . Its fermionic component is given in terms of radicals and trigonometric functions. 
The classical modified equation (11.211) possesses a solution (I2.15P expressed in terms of 
elliptic functions, whose imaginary part can be used to approximate a bump function. 
This solution also appears as the bosonic component of a solution (I5.10p and 05. 11 j) of 
the supersymmetric equations (14.81) and (14. 9 p with e — — 1. For the invariant solutions 
of both supersymmetric extensions, there are instances (eg. ( 15.181) . ( 15.191) and ( 15.251) ) 
where the fermionic field ip is an arbitrary function of a single quantity. This permits 
us consider a wide range of possibilities for the fermionic medium. Equations (13.71) and 
(13. 8p represent the static case which involves two spatial variables. It is expected that a 
supersymmetric extension of the Gaussian fluid flow equations in (2+ 1) dimensions may 
yield richer classes of physically relevant solutions. 
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